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Background: BRAF is mutated in -42% of human melanomas (COSMIC. http://www.sanger.ac.uk/genetics/CGP/ 
cosmic/) and pharmacological BRAF inhibitors such as vemurafenib and dabrafenib achieve dramatic responses in 
patients whose tumours harbour BRAF V600 mutations. Objective responses occur in -50% of patients and disease stabil- 
isation in a further -30%, but -20% of patients present primary or innate resistance and do not respond. Here, we investi- 
gated the underlying cause of treatment failure in a patient with BRAF mutant melanoma who presented primary 
resistance. 

Methods: We carried out whole-genome sequencing and single nucleotide polymorphism (SNP) array analysis of five 
metastatic tumours from the patient. We validated mechanisms of resistance in a cell line derived from the patient's 
tumour. 

Results: We observed that the majority of the single-nucleotide variants identified were shared across all tumour sites, 
but also saw site-specific copy-number alterations in discrete cell populations at different sites. We found that two ubiqui- 
tous mutations mediated resistance to BRAF inhibition in these tumours. A mutation in GNAQ sustained mitogen-acti- 
vated protein kinase (MAPK) signalling, whereas a mutation in PTEN activated the PI3 K/AKT pathway. Inhibition of both 
pathways synergised to block the growth of the cells. 

Conclusions: Our analyses show that the five metastases arose from a common progenitor and acquired additional 
alterations after disease dissemination. We demonstrate that a distinct combination of mutations mediated primary re- 
sistance to BRAF inhibition in this patient. These mutations were present in all five tumours and in a tumour sample 
taken before BRAF inhibitor treatment was administered. Inhibition of both pathways was required to block tumour cell 
growth, suggesting that combined targeting of these pathways could have been a valid therapeutic approach for this 
patient. 
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introduction 

Melanoma is a potentially deadly form of skin cancer. The 
metastatic form of the disease has a poor prognosis with a 
median survival of 6-9 months and an overall survival of 10%- 
15% [1]. BRAF, which encodes a protein kinase, is mutated in 
-42% of human melanomas. Drugs that inhibit BRAF, such as 
vemurafenib and dabrafenib, can achieve dramatic responses in 
patients whose tumours harbour V600 BRAF mutations, with 
objective responses in ~50% of patients and disease stabilisation 
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in a further ~30% [2, 3]. Thus, BRAF is a validated therapeutic 
target in melanoma, but despite these impressive results, 
acquired (secondary) resistance is almost universal, even when 
the initial response leads to profound tumour regression. In 
most cases, acquired resistance to BRAF inhibitors is due to 
MAPK pathway reactivation and is driven by events that can 
occur upstream, downstream or even at the level of BRAF [4, 5]. 
In the remainder of cases, resistance is driven by mitogen- 
activated protein kinase /extracellular signal-regulated kinase 
(MEK/ERK) -independent mechanisms. Notably, inter- and 
intra-lesional molecular heterogeneity can also shape acquired re- 
sistance to BRAF inhibitors, resulting in different mechanisms of 
resistance in individual tumours from a single patient [6, 7] or 
the co- existence of sensitive and resistant clones within a single 
tumour [8] . 
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In addition to treatment failures due to acquired resistance, 
~20% of patients present intrinsic (primary/innate) resistance 
and do not derive any benefit from BRAF inhibitors [2, 3]. 

However, while impressive advances have been made in our 
knowledge of mechanisms of acquired resistance, there is limited 
mechanistic insight into intrinsic resistance. RAC1 P29S mutations 
have recently been reported in one patient with intrinsic resist- 
ance and two patients with transiently stable disease followed by 
disease progression after administration of a BRAF inhibitor [9] 
but the functional significance of these intriguing observations 
was not elucidated. We used whole-genome sequencing (WGS) 
to characterise five metastatic tumours from a BRAF mutant mel- 
anoma patient who presented intrinsic resistance. We observed 
evidence of inter-lesional heterogeneity that pointed to clonal 
complexity, but identified concomitant mutations in GNAQ and 
PTEN as the mechanisms of resistance in all five tumours. Thus, 
we describe complex but ubiquitous mechanisms of intrinsic 
resistance that explain treatment failure in this patient. 

methods 

tissue and blood collection 

Following patient consent, tumour tissue and blood were obtained at surgery 
and in the outpatients department, respectively. All samples were anon- 
ymised and access to samples and clinical data was restricted in accordance 
with the Human Tissue Act and Multi- centre Research Ethics Committees 
(MREC) guidelines. All the studies were conducted in accordance with a 
study protocol CCR3097, approved by the Royal Marsden Hospital Research 
Ethics Committee on 21 October 2008. 

DNA extraction 

The tissue was disrupted using Precellys®24 tissue homogeniser (Precellys, 
Dublin, Ireland). Subsequently, the tissue was lysed in ATL buffer (Qiagen, 
Manchester, UK) and Proteinase K (Qiagen) for 24 h at 56°C. The DNA was 
purified using the DNAaesy Blood and Tissue spin-column Kit (Qiagen) 
according to manufacturer's instructions. DNA from bufiy coat and cultured 
cells was isolated using the same kit. DNA was quantified using the 
PicoGreen dsDNA Quantification Reagent (Invitrogen) according to manu- 
facturer's recommendations. The structural integrity of the tumour DNA 
was assessed by gel electrophoresis. 

Sanger sequencing 

Regions of interest were amplified by PCR. The products were directly 
sequenced using dye-terminator chemistry. Samples were analysed on a 
3100 Genetic Analyser (Applied Biosystems, Paisley, UK). Sequences were 
visualised using Sequencher software (Gene Codes Corporation, Ann Arbor, 
MI, USA). 

WGS and somatic variant detection 

Sequencing of the matched normal and tumour samples was carried out 
using unchained combinatorial probe anchor ligation chemistry on arrays of 
self- assembling DNA nanoballs [10]. The reads were aligned to the NCBI 
build 37 reference genome. The gross mapping yield, haploid coverage and 
percentage of fully called bases for each genome are reported in supplemen- 
tary Table SI, available at Annals of Oncology online. Variations between the 
reference genome (NCBI build 37) and each of the samples were called and 
scored using a local de novo assembly algorithm [10]. 

Somatic single-nucleotide variants (SNVs) and indels between the normal 
and tumour samples uncovered using the calldiff function of cgatools (http:// 
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cgatools.sourceforge.net/), which assigns a somatic score to each SNV and 
additionally a somatic score to short insertions and deletions. The somatic 
score is a measure of the confidence that each SNV/indel is a true somatic 
variant. Somatic variants with a somatic score >— 10 were considered high- 
confidence variants. Variants present in single nucleotide polymorphism 
database (dbSNP) were excluded from further analysis. SNVs with somatic 
scores <-10 in a genome were designated true somatic variants if present in 
one of the other genomes with a somatic score >-10. A somatic score 
threshold of >0 was used to identify high- confidence site-specific SNVs. 

Somatic structural variations (SVs) in the primary and metastatic tumours 
were called using the junctiondiff function of cgatools. High -confidence SV 
junctions were those that had at least 10 mate-pairs in a cluster, in which de 
novo assembly of the junction was successful, had a high mapping diversity 
and for which there was an absence of specific repeat sequences on the left and 
right side of the junction. 

copy number alterations 

SNP array. Genome wide genotyping was conducted on the Affymetrix 
Genome- Wide Human SNP Array 6.0. Genotypes were called using the apt- 
probes et- genotype algorithm (birdseed-v2 method) of Affymetrix Power 
Tools. logR ratios (LRRs) and B-allele frequencies (BAFs) were estimated by 
PennCNV-Afty and allele- specfic copy-number alterations between the 
germline and tumour genomes were called by the ASCAT algorithm using 
default parameters [11]. 

whole-genome sequencing. Read depth was corrected for GC content 
and normalised for average haploid genome coverage. Log2 ratios of the tumour 
to normal average read depths for windows of 2 kb were calculated. These values 
were smoothed and segmented as implemented in the BioConductor package 
DNAcopy using default parameters [12]. Bootstrapped, unsupervised clustering 
of the WGS CNA data was conducted using the R package pvclust (http://www. 
is.titech.ac.jp/ -shimo/ prog/ pvclust/) . 

estimation of tumour purity 

We conducted an analysis of the purity of the tumours using the SNP array 
data with the ASCAT algorithm [11]. This algorithm reports tumour cell frac- 
tions of 69%, 62%, 59%, 80% and 66% for sites 1, 2, 3, 4 and 5, respectively. 

development of a continuous cell line 
from patient's tumour 

Tumour tissue was minced using sterile crossed blades under sterile condi- 
tions. The pieces were plated sparsely in 6-well uncoated or collagen coated 
plates. Enough culture media was added to cover each tissue piece and the 
plates incubated in a humidified atmosphere at 37°C and 5% C0 2 and either 
20% 0 2 or 3% 0 2 . When cells emerging from tissue pieces became confluent, 
they were passaged at dilution of 1:2. 

immunofluorescence of cultured cells 

Cells were grown on cover slips in 24-well plates and fixed with 4% formal- 
dehyde in phosphate buffer solution (PBS) for 15 min. The cells rinsed in 
PBS three times, for 5 min each. Blocking buffer (1% bovine serum albumin 
(BSA) in PBS) was applied for 60 min, following which the specimen was 
incubated with the primary antibody overnight at 4°C. The specimen was 
rinsed in PBS three times for 5 min each, and the secondary antibody 
applied for 2 h at room temperature. The secondary antibody was rinsed off 
with PBS once and DAPI once, followed by three PBS washes, 5 min each. 
Coverslips were mounted with Vectashield and slides sealed with clear nail 
varnish. The results were visualised in a confocal microscope. 
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reagents 

For western blotting, the following antibodies were used: rabbit anti-ERK2 
(C-14) (Santa Cruz Biotechnology, Heidelberg, Germany); mouse anti- 
Tubulin, and mouse anti-ppERKl/2 (Sigma, Dorset, UK); rabbit anti- 
PTEN and rabbit anti-ppAKT (Cell Signalling, Hitchin, UK). For 
immunofluorescence, the following antibodies were used: rabbit anti-SlOO 
(Dako, Ely, UK) and mouse anti-HMB45 (Abeam, Cambridge, UK). 
MK2206 was obtained from Stratech Laboratories (Suffolk, UK). 
PD 1843 52 and PLX4720 were synthesised in-house. All drugs were pre- 
pared in Dimethyl sulfoxide (DMSO). Fluorescence in situ hybridisation 
(FISH) probes for chromosome 3, 7 and 17 were obtained from UroVysion, 
Abbott Molecular. 

cell culture techniques 

A375P cells were cultured in DMEM supplemented with 10% foetal bovine 
serum. A total of 22 092 cells were cultured in MCDB153 supplemented 
with 20% Leibovitz's L-15, 2% foetal bovine serum, insulin (5 (ig/ml) and 
CaCl2 (1.68 mM). A375P cells were incubated at 37°C, 10% C0 2 and 20% 
0 2 . A total of 22 092 cells were incubated at 37°C, 5% C0 2 and 3% 0 2 . For 
inhibitor treatment, the drugs were dissolved in DMSO and added to the 
medium for 3 h, unless otherwise stated. 

cell proliferation assays 

96-well plate format. Cells were seeded at a density of 1000-5000 cells 
per well in 100 ul medium. Cells were treated with drug in triplicate for 3-10 
days. For the longer incubations (>5 days), the media and the drugs were 
replenished every 3 days. Cell growth was measured using the CellTiter-Glo 
Luminescent Cell Viability Assay (Promega, Southampton, UK). The 
luminescence was recorded on a luminometer (SpectraMax M5, Molecular 
Devices, Wokingham, UK). 

6-well plate format. Cells were seeded at a density oflO 000-20 OOOcells 
per well in 1000 ul medium. Cells were treated with drug for 10 days. The 
media and drug were replenished every 3 days. The media were removed and 
cells rinsed with PBS. Two millilitres of a mixture of 6% glutaraldehyde and 
0.5% crystal violet were added to each well and incubated for 30 min. Crystal 
violet staining was quantified using ImageJ software (www.rsb.info.nih.gov/ij). 

Spheroids assays. One hundred microlitre of 50% Matrigel (diluted with 
medium) was placed into each well of a 96-well plate. The plate was incubated 
at 37°C for 1 h to set the matrigel and cells (1000/well) were seeded on top of 
the matrigel plug. Drugs were added at the time of cell seeding and both media 
and drugs were replenished every 2 days. Ten to 14 days after seeding images 
were taken for quantification of spheroid size and branching. At least four 
representative images were taken of each well, and ImageJ software was used 
for the analysis. 

small interfering RNA transfection. 2.5 x 10 5 cells were seeded per 
well in a 6-well dish and incubated overnight. The following day the cells 
were transfected with Lipofectamine (Invitrogen, Paisley, UK) and 20 nM of 
the following small interfering RNA (siRNA) pools: sc-35429 (Santa-Cruz), 
s5886 (Life Technologies, Paisley, UK) and s5887 (Life Technologies), or 
scrambled control. Cells were incubated for 48 h before preparation of cell 
lysates, or re-plating for a proliferations assay. 

expression constructs 

A plasmid with the entire GNAQ coding region was purchased from Open 
Biosystems. The Q209P mutation was created by standard PCR-directed mu- 
tagenesis approaches. A Spel/SphI restriction fragment was released and 



introduced into a similarly digested pMCEF- vector. pMCEF- was derived 
from pMClneo [13]. 

DNA transfections 

For transient transfection of A375P cells, 2.5 x 10 5 cells were seeded per well 
in a 6-well dish and incubated overnight. Lipofectamine (Invitrogen) was 
used according to manufacturer's instructions. Cell lysates were prepared 2 
days following transfection. For generation of stable cell lines, A375P cells 
were transfected with the various vectors carrying a NeoR cassette using 
Lipofectamine and selected with G418 (Invitrogen). The cells were harvested 
for isolation of mRNA and positive clones selected by quantitative real-time 
PCR (qRT-PCR). 

quantitative real-time PCR 

Primers for each candidate gene were designed off 3' regions of the cDNA 
sequences. GNAQ probe (Hs00387073_ml) and the GAPDH Endogenous 
Control probe were obtained from Applied Biosystems. For each reaction, 
10 ul of Precision Mastermix (PrimerDesign), 1 ul of the probe and 25 ng of 
the cDNA template were added and made up with to 20 ul nuclease-free 
H 2 0. All samples were tested in triplicates and H 2 0 was used as a negative 
control. qRT-PCR was carried out on an Applied Biosystems 7900HT Fast 
Real Time Machine (Applied Biosystems). Relative expression was calculated 
using the AAQ method with GAPDH as an internal control. 

data access 

WGS and SNP array data will be available from the European Genome- 
phenome Archive under study accession number EGAS0000100058. 

results 

A 26-year-old male presented with a thick ulcerated cutaneous 
melanoma that, despite surgical resection, metastasised to the 
mediastinal lymph nodes, so systemic therapy was considered. 
Genotyping of archival tumour tissue revealed a T>A, p.V600E 
mutation in BRAF, and the patient was enrolled into a clinical 
trial with vemurafenib (BRIM-3) [14]. However, radiological as- 
sessment at 8 weeks revealed disease progression with a fivefold 
increase in the volume of the mediastinal mass and so vemurafe- 
nib was discontinued (Figure 1A). Five weeks later, a new sub- 
cutaneous tumour appeared and was removed (site 1; Figure 1A 
and B), and fresh tissue obtained for study. Two cycles of decar- 
bazine and cisplatin were administered, but the patient derived 
no benefit. Seven weeks after completion of chemotherapy a 
further four subcutaneous tumours were removed (sites 2, 3, 4, 5; 
Figure 1 A and B) and we obtained fresh tissue for study. A forma- 
lin-fixed subcutaneous metastasis removed 1 month before 
vemurafenib therapy was available as a baseline sample. 

We confirmed that the T>A, p.V600E mutation in BRAF was 
present in all five tumours (supplementary Figure SI, available 
at Annals of Oncology online), so to investigate mechanisms of 
resistance to vemurafenib, we carried out WGS using the 
Complete Genomics platform [10] (supplementary Table SI, 
available at Annals of Oncology online). Compared with the 
patient's germline DNA, we observed ~40 000 predicted 
somatic SNVs per tumour (supplementary Table SI, available at 
Annals of Oncology online), a level that is comparable with the 
mutation burden in other sun-exposed melanomas [15, 16]. The 
majority (36 507; 77%) of the SNVs were present in all tumours 
(Figure 2 A) and decarbazine/cisplatin treatment did not 
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Figure 1. Clinical case under study. (A) The clinical course of the patient is depicted over a period of 12 months, with time points and duration of different 
therapeutic interventions indicated. Metastases from sites T, '2', '3', '4' and '5' were obtained as fresh tissue. (B) Diagram showing the site of the subcutaneous 
sites from which metastatic tissue was removed. The tumour from site 1 was propagated as a continuous cell line. 



significantly increase SNV load in tumours #2, #3, #4 and #5 
(supplementary Table SI, available at Annals of Oncology 
online). Approximately 85% of the common somatic SNVs were 
C>T (G>A) transitions (Figure 2B; supplementary Table SI, 
available at Annals of Oncology online), of which over 90% were 
at the 3' end of pyrimidine dinucleotides (Figure 2C, 'All sites'; 
supplementary Table SI, available at Annals of Oncology 
online), consistent with a UV- induced DNA damage signature 
[17]. Notably, the high- confidence private SNVs from the indi- 
vidual tumours did not bear this signature (Figure 2C, sites 1-5; 
Figure 2D). Our data identified 283 predicted non-synonymous 
(ns) SNVs in one or more tumours, 219 of which were common 
to all five sites. Sanger sequencing revealed that the majority of 
the non-ubiquitous nsSNVs were false positives or miscalled 



germline variants (supplementary Table S2, available at Annals 
of Oncology online), and no nsSNVs were confirmed as exclu- 
sive. Genome wide, we predicted ~424 somatic small deletions 
and ~430 somatic small insertions, only 295 of which were 
present in all tumours. Sanger sequencing of these mutations 
confirmed only two coding region indels, both of which were 
present in all five tumours (supplementary Table S3, available at 
Annals of Oncology online). 

We used WGS and SNP array analysis [11] to reveal somatic 
SVs and copy-number alterations (CNAs). The SV analysis pre- 
dicted 89-129 high- confidence structural variants, 29 of which 
were common to all tumours (supplementary Figure S2 and 
Table S4, available at Annals of Oncology online). The CNA ana- 
lysis revealed very similar landscapes for sites 1, 2 and 3, 
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Figure 2. Genomic analysis of tumours. (A) Venn diagram showing numerical overlap in the high- confidence SNVs in tumour sites 1-5. (B) Graph showing 
mutation signature of the 36 507 common SNVs. (C) Stacked histograms showing the mutation context of the common C to T (G to A) transitions in all 
five tumours (All sites; top left panel), and for the high- confidence (somatic score >0) private C to T (G to A) transitions in each individual tumour (sites 1-5). 
The proportion of each nucleotide ±1 bp of the mutation is shown (A: orange, C: pink, G: blue, T: green). (D) Graph showing mutation signature of the high- 
confidence (somatic score >0) private (site-specific) SNVs for each individual tumour (sites 1-5). (E) Genome plots of somatic copy-number alterations in the site 
1 (top panel) and site 5 (bottom panel) tumours as determined by normalised coverage across 2 kb windows. (F) Figure showing bootstrapped unsupervised clus- 
tering of the WGS copy-number alteration data. The approximately unbiased (au; left) P- value and bootstrap probability (bp; right) values are indicated. 



whereas sites 4 and 5 showed higher levels of aneuploidy 
(Figure 2E; supplementary Figure S3, available at Annals of 
Oncology online; supplementary Table S5, available at Annals of 
Oncology online). In particular, site 5 showed increased copy 
number of chromosome 7 (supplementary Figure S3, available 
at Annals of Oncology online), and this tumour presented higher 
levels of nuclear polymorphism (supplementary Figure S4, 
available at Annals of Oncology online, supplementary Table S6, 
available at Annals of Oncology online). Bootstrapped unsuper- 
vised clustering of the WGS CNA data revealed a phylogeny in 
which sites 1, 2 and 3 segregated from sites 4 and 5 (Figure 2F). 
We carried out FISH, and this revealed that the cells in sites 2 
and 3 were homogeneous and had largely diploid nuclei, 
whereas those in sites 4 and 5 were aneuploid (supplementary 
Table S5, available at Annals of Oncology online). Specifically, 
site 4 presented three to five chromosomes per nucleus, and site 
5 was comprised of a mixed cell population of small cells with 
predominantly diploid nuclei and large cells with up to 15 chro- 
mosomes per nucleus (supplementary Table S6, available at 
Annals of Oncology online). Thus, all tumours shared several 
common events, including being hemizygous for chromosome 
10 (Figure 2E; supplementary Figure S3, available at Annals of 
Oncology online), but more recently, the heterogeneous cells in 



sites 4 and 5 appeared to have acquired additional chromosomal 
gains. 

We established a continuous cell line (designated 22 092 cells) 
from tumour #1 and confirmed that the T>A, p.V600E muta- 
tion in BRAF was present in this line (supplementary 
Figure S5A, available at Annals of Oncology online). The cells 
stained positive for the melanoma markers HMB45 and SI 00 
(supplementary Figure S5B, available at Annals of Oncology 
online), and the majority of the SNVs identified in tumour #1 
persisted in this line (supplementary Figure S5C, available at 
Annals of Oncology online). Notably, compared with A375 mel- 
anoma cells, which also express BRAF V600E , the 22 092 cells 
were significantly less sensitive to the BRAF inhibitor PLX4720, 
but displayed equipotency to the MEK inhibitor PD 1843 52 
(Figure 3 A). Accordingly, ERK phosphorylation was inhibited 
by PLX4720 in A375 cells, but persisted in 22 092 cells, whereas 
PD184352 inhibited ERK equally in both lines (Figure 3B). 
Thus, despite the BRAF V600E mutation, 22 092 cells were resist- 
ant to BRAF inhibitors. 

No mutations were observed in RAS, MEK or ERK, but all five 
tumours harboured an A>C, p.Q209P mutation in GNAQ (sup- 
plementary Figure S6, available at Annals of Oncology online). 
Sanger sequencing of the pre-treatment lesion confirmed that this 
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Figure 3. Oncogenic GNAQ contributes to BRAFi resistance. (A) Graph showing A375P and 22 092 cell growth in the presence of increasing concentrations 
of PLX4720 and PD184352. (B) Western blots for phospho-ERK (ppERK) and total ERK2 (loading control) in A375P and 22 092 cells treated with the indi- 
cated concentrations in micromolar of PLX4720 or PD184352 for 3 h. (C) Western blot for phospho-ERK (ppERK) and tubulin (loading control) in A375P 
cells expressing wild-type (GNAQ WT ) or mutant (GNAQ Q209P ) GNAQ and treated with the indicated concentrations of PLX4720 for 3 h. (D) Graph showing 
growth of A375P cells expressing wild-type (GNAQ WT ) or mutant (GNAQ Q209P ) GNAQ and treated with 1 uM PLX4720 for 10 days. Cell growth was quan- 
tified by crystal violet. (E) Western blot for phospho-ERK (ppERK) and tubulin (loading control) in A375P and 22 092 cells transfected with scrambled control 
(Control), or three different GNAQ siRNAs (siRNA 1, siRNA 2, siRNA 3) targeting GNAQ for 48 h. ppERK levels were quantified using ImageJ and the mean 
reduction (± SD) in levels is shown. (F) Graph showing growth of 22092 cells transfected with scrambled control (siRNA scr) or two GNAQ siRNAs (siRNA 1, 
siRNA 2) for 48 h and then treated with PLX4720 for 72 h. (G) Graph showing growth of A375P cells transfected with scrambled control (siRNA scr) or two 
GNAQ (siRNA 1, siRNA 2) siRNAs for 48 h and then treated with PLX4720 for 72 h. 
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mutation was present before vemurafenib therapy (supplemen- 
tary Figure S6, available at Annals of Oncology online). Activating 
mutations in GNAQ occur in ~50% of uveal melanomas [18] but 
are rare in cutaneous melanoma [18, 19]. To date, three cases of 
concurrent mutations in GNAQ and BRAF have been reported 
in melanoma (http://www.cbioportal.org/), but the mutations in 
GNAQ (D333fs, G64A and P185S/A302G) were not at the 
common Q209 hotspot, so the functional significance of those 
alterations is unclear. However, oncogenic GNAQ activates ERK 
in uveal melanoma [18, 19, 20], so we investigated its effects in 
BRAF-mutant melanoma cells. When GNAQ Q209P was expressed 
in A375 cells (supplementary Figure S7, available at Annals of 
Oncology online), it increased basal ERK phosphorylation and 
reduced ERK sensitivity to PLX4720 (Figure 3C) and, critically, it 
also reduced the sensitivity of A375 cells to PLX4720 (Figure 3D). 
Conversely, when GNAQ Q209P was depleted in 22 092 cells using 
siRNA (supplementary Figure S8, available at Annals of Oncology 
online), ERK activity was suppressed, whereas depletion of wild- 
type GNAQ from A3 75 cells had less effect on ERK activity 
(Figure 3E). Notably, depletion of GNAQ Q209P increased the sen- 
sitivity of 22 092 cells to PLX4720 (Figure 3F), whereas depletion 
of wild-type GNAQ did not increase A375 cell sensitivity to this 
drug (Figure 3G). 

Intriguingly, one of the other confirmed mutations in all of 
the tumour samples, including the pre-treatment sample, was a 
4-bp deletion in PTEN c.950_953delTACT (supplementary 
Figure S9, available at Annals of Oncology online) that introduces 
a frameshift at codon V317 and a stop at codon 319 within the 
C2 binding domain. Sanger sequencing revealed wild-type PTEN 
sequence traces in the tumours, but not the 22 092 cells (supple- 
mentary Figure S8, available at Annals of Oncology online), sug- 
gesting that the wild-type PTEN sequence in the tumour samples 
is the result of stromal contamination. Accordingly, we observed 
chromosome 10 loss in all tumours (Figure 2E; supplementary 
Figure S3, available at Annals of Oncology online) and observed 
that, in the 22 092 cells, PTEN was not expressed and AKT phos- 
phorylation was elevated (Figure 4A). The allosteric AKT inhibi- 
tor MK2206 inhibited AKT phosphorylation in 22 092 cells 
(Figure 4B) and 22 092 cells were more sensitive to MK2206 than 
A3 75 cells (Figure 4C). Taken together, our data show that, in the 
tumours from this patient, ERK activity was elevated due to the 
GNAQ P209 mutation, and AKT activity was elevated by loss of 
PTEN expression. Accordingly, PD 184352 and MK2206 syner- 
gised to inhibit the growth of 22 092 cells in long-term prolifer- 
ation assays (Figure 4D), and when the cells were grown in 
suspension (Figure 4E). 

discussion 

Here, we analysed the genomes of five metastatic BRAF V600E 
melanomas from a patient who presented intrinsic resistance to 
vemurafenib. Previous studies have reported discordance in 
BRAF mutation status in different tumours from other patients 
[21] and this could drive differential response to BRAF inhibi- 
tors. However, we confirmed that the T>A, p.V600E mutation 
in BRAF persisted in all five sites and that it was present in the 
pre-treatment sample. Our WGS revealed an A>C, p.Q209P 
mutation in GNAQ that was also present in all five tumours. We 
demonstrated that GNAQ Q209P sustained ERK activity in BRAF 
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mutant melanoma cells in the presence of a BRAF inhibitor, 
allowing the cells to grow even when BRAF V600E was inhibited. 
Our sequencing also revealed a ubiquitous PTEN frame-shift de- 
letion that was associated with AKT hyper- activation. Both 
GNAQ and PTEN mutations were detected in the pre-vemura- 
fenib sample, observations that are consistent with features of 
intrinsic mechanism of resistance; however, since the first as- 
sessment was carried out at 8 weeks after starting treatment, we 
cannot exclude the possibility that the patient had a transient 
early response followed by rapid progression. 

It has been reported that mutations in the RAS small G-proteins 
can mediate acquired resistance to BRAF inhibitors in BRAF 
mutant melanoma [6] and that GNAQ-mutant uveal melanomas 
do not respond to BRAF inhibitors [22]. However, mutations in 
GNAQ have not previously been reported to mediate resistance to 
BRAF inhibitors in BRAF-mutant melanoma. BRAF and GNAQ 
mutations are thought to be mutually exclusive and activating 
GNAQ mutations have not been reported to occur concomitant 
with BRAF mutations. It has been reported previously that PTEN 
mutations can mediate acquired resistance to BRAF inhibition in 
melanoma pre-clinical models [23, 24]. Although down-regulation 
of PTEN has an association with overall lower responses and 
shorter progression-free interval in patients [25, 26], the presence 
of PTEN alterations does not always preclude a response to BRAF 
inhibitors [9]. Critically, in our study, we show that both the 
MEK/ERK and AKT pathways needed to be targeted to mediate 
effective inhibition of proliferation of cells derived from this 
patient's tumour, suggesting that loss of PTEN did contribute to 
intrinsic resistance to BRAF inhibition in this patient. 

An important consequence of the clonal nature of cancers is 
spatial and temporal variation in tumour composition [27]. 
Spatial heterogeneity has been demonstrated in several tumour 
types, including renal [28] and pancreatic [29] cancers. We 
observed a high degree of homogeneity at the single-nucleotide 
level in the tumours from this patient, but nevertheless observed 
more CNAs in the tumours at sites 4 and 5. Furthermore, at site 5 
we observed a mixture of cells with distinct genomic and pheno- 
typic features. This suggests that, while the tumours all evolved 
from a common metastatic progenitor, the cancer cells at sites 4 
and 5 acquired additional genomic alterations after disease dis- 
semination. We note that all of the tumours had lost a copy of 
chromosome 18 and that loss of 18q is linked to chromosomal in- 
stability [30]. Furthermore, although the DNA damaging agents 
cisplatin and decarbazine did not increase SNV burden or alter 
the mutation signature in sites 2, 3, 4 and 5 compared with site 1, 
we cannot exclude the possibility that these agents caused the 
additional chromosomal alterations we observed at sites 4 and 
5. Critically, these findings of intra- and inter-lesional heterogen- 
eity caution against the sampling bias that can occur when 
genomic data from single- tumour biopsies are analysed. 

Our data reveal that resistance in this patient was not 
mediated by acquired mutations or sub-clonal evolution under 
selection pressure exerted by treatment, but by a combination of 
pre-existing mutations in two genes that activate at least two 
pro-growth/pro-survival pathways. These mutations were 
present in all tumours and predated vemurafenib administra- 
tion, explaining why the patient presented intrinsic resistance. 
Two recent studies have also reported multiple mechanisms of 
acquired resistance to BRAF inhibitors within the same tumour 
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Figure 4. Loss of PTEN contributes to BRAFi resistance in patient cells (A) Western blot for phospho-AKT (ppAKT), PTEN and total ERK2 (loading 
control) in A375P and 22 092 cells. (B) Western blot for ppAKT, phospho-ERK (ppERK) and total ERK2 (loading control) in 22 092 cells treated with the indi- 
cated concentrations of MK2206 for 3 h. (C) Graph showing growth of A375P and 22 092 treated with increasing concentrations of MK2206 for 10 days. 
(D) Photograph showing crystal violet staining of 22 092 cells treated with DMSO, MK2206 (1 (iM), PD184352 (0.1 (iM), or both for 10 days. (E) Graph 
showing growth of 22 092 spheroids in collagen following treatment with the indicated concentrations of PD 184352, MK2206 or both for 7 days. 
Representative images of the treated spheroids are shown below and were assessed by measuring individual spheroids with ImageJ software in triplicate for each 
condition. **P value <0.01 individual treatments compared with DMSO -treated control; ***P value <0.001 combination treatment compared with PD 184352 
alone. 
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biopsy [9, 31]. Together with the findings reported here, these 
observations highlight the genetic complexity that can drive 
drug resistance in melanoma patients. 

We demonstrated that AKT and MEK inhibitors synergised 
to block the growth of the tumour cells, suggesting that com- 
bined inhibition of these pathways may have been an effective 
approach for this patient. Our findings show how comprehen- 
sive tumour profiling combined with knowledge of signalling 
biology can identify biomarkers that may predict treatment re- 
sistance and may also reveal effective drug combinations. 
Clearly, this will help formulate more effective therapeutic strat- 
egies for individual patients. 
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